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ABSTRACT 

It is difficult to imagine a planet formation model that does not at some 
stage include a gravitationally unstable disc. Initially unstable gas-dust discs 
may form planets directly, but the high surface density required has motivated 
the alternative that gravitational instability occurs in a dust sub-layer only after 
grains have grown large enough by electrostatic sticking. Although such growth 
up to the instability stage is efficient for laminar discs, concern has mounted as 
to whether realistic disc turbulence catastrophically increases the settling time, 
thereby requiring additional processes to facilitate planet formation on the needed 
time scales. To evaluate this concern, we develop a model for grain growth that 
accounts for the influence of turbulence on the collisional velocity of grains and on 
the scale height of the dust layer. The relative effect on these quantities depends 
on the grain size. The model produces a disc-radius dependent time scale to 
reach a gravitationally unstable phase of planet formation. For a range of dust 
sticking and disc parameters, we find that for viscosity parameters a < 10~ 3 , 
this time scale is short enough over a significant range in radii R that turbulence 
does not catastrophically slow the early phases of planet formation, even in the 
absence of agglomeration enhancement agents like vortices. 

Subject headings: accretion discs - planetary systems: formation - planetary 
systems: protoplanetary discs 



1. Introduction 

Planets are believed to form in the gas and dust discs that surround newly formed 
stars. In principle, planets can form purely from gravity, via direct gravitational instability 
in the initial circumstellar disc (Boss 1997). However, this mechanism of planet formation 
requires higher initial densities than commonly presumed. For a standard MMSN (minimum 
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mass solar nebula) model for the gas and dust discs around a young star, both the initial 
dust and gas discs are gravitationally stable and the only force available for the early proto- 
planetesimal stage of planet formation is electrostatic sticking. For planets to form, the dust 
must therefore grow from electrostatic forces until gravity can play a significant role. 

The basic model for how such a state could can arise was put forth by Goldreich and 
Ward (Goldreich & Ward (1973), hereafter GW). In their proposed route to planet formation, 
dust grains initially collide and stick. The growing grain mass eventually reduces the dust's 
thermal velocity dispersion and therefore, the scale height of the dust disc. Eventually the 
grains settle into the disc's mid-plane until the critical density at which this dust disc becomes 
gravitationally unstable is reached and the formation of kilometer sized planetesimals is 
facilitated. Subsequent accretion of surrounding gas would then complete the core-accretion 
model for giant planet formation. Recent observations of CoKu Tau/4 (D'Alessio et al. 
2005) seem to imply that massive planets must be able to form within ~ 10 6 yr, which is the 
tightest constraint to date on the total available time for any planet formation mechanism. 

In the absence of turbulence in the gas disc, the GW model is very efficient and the time 
scale to grow grains to the size at which enough settling occurs for gravitational instability 
to occur is a small fraction of 10 6 years. However, Weidenschilling (Weidenschilling 1980) 
argued that turbulence can catastrophically prevent the required early growth phase by stir- 
ring up the dust disc, thereby delaying or preventing the subsequent gravitational instability. 
Because of prevalent sources of disc turbulence (such as dust-gas shear, Cuzzi et al. (1993), 
Champney et al. (1995) and MRI instability Balbus & Hawley (1991)) a turbulent disc is 
likely the norm rather than the exception and the GW theory must be revisited. Ironically, 
even turbulence driven from the dust settling itself (Ishitsu & Sekiya 2003), might prema- 
turely quench the GW process. The potential show-stopping effect of turbulence has led to 
a substantial body of work incorporating additional physics, such as anti-cyclonic vortices, 
that can accelerate the agglomeration of dust grains should the early GW phase fail. But 
even if extra processes are present and helpful, the need for such processes has remained 
unclear. Despite the conceptual concerns induced by turbulence, its actual effect on planet 
growth has not been conclusively calculated. Regardless of the details of the turbulence, 
there will exist a dust grain size for which the effects of the turbulence on the dust are 
sufficiently weak that the GW process can proceed to instability. The question therefore is 
whether the presence of turbulence catastrophically slows or even stops grain growth before 
that size can be reached. We attempt to answer that question. 

The effect of turbulence on dust growth is two-fold. On the one hand it can increase 
the dust's collisional velocities (and so the collisional rate) which is a positive effect for 
growth. Some early work on this was done by Voelk et al. (1980) (hereafter VMRJ), further 
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developed by Markiewicz et al. (1991) (hereafter MMV) and Cuzzi & Hogan (2003). On the 
other hand, it increases the scale height of the dust (Dubrulle et al. (1995), who used the 
results of VMRJ) which acts both to decrease the density (and so the collisional rate) and 
to prevent instability. Understanding the net effect of turbulence amounts to understanding 
its relative effect on these two parameters as a function of grain size and radial location. 
Suchs efforts have been made by Kornet et al. (2001) and Stepinski & Valageas (1997) using 
models by Cuzzi et al. (1993). In this paper we develop a different model than those of 
VMRJ or Cuzzi et al. (1993) for the effects of turbulence on grain growth and use it to 
calculate the time evolution of grain growth as a function of dust size, disc location and 
turbulent strength; so as to determine how long it takes for the grains to grow large enough 
to settle in the disc and reach the gravitationally unstable density. Both our approach to 
calculating the effects of turbulence and our results differ from the previous work as will be 
discussed in section 3.3. We do not discuss the physics of the gravitationally unstable and 
post-gravitationally unstable regimes. 

In section 2 we derive the basic equation for grain growth, describe the disc model and 
show how turbulence and drag on grain motion are included. In terms of the parameters 
we develop, the grain growth equation depends on the (dimensional) ratio of dust collisional 
velocity and the dust scale height. In section 3 we derive explicit formulae for these key 
parameters. In section 4 we show that there are 5 regimes for the combination of these 
parameters and discuss how grain growth proceeds as a progression through these regimes. 
In section 5 we collect key results for the final grain sizes and time scales needed to get to 
the scales at which the disc becomes gravitationally unstable for fully turbulent discs and for 
discs with dead zones. In section 6 we discuss the implications of our results for constraining 
the strength of turbulence for which rapid planet growth can still occur. In section 7 we 
examine our key assumptions and we conclude in section 8. 



2. Basic model of dust grain growth in a turbulent disc 

We consider a protoplanetary disc to be a dust disc embedded in a turbulent gas disc. 
At any one time, we approximate all dust grains as spherical with identical radii r^, constant 
density pd and hence dust-dust collisional cross-section cr<ta — ^ nr d an d mass rrid = |7rp^r^. 
Although the dust in a real disc would have a spectrum of sizes, the Central Limit theorem 
implies that, in the absence of significant destruction from grain-grain collisions, the dust 
spectrum will rapidly approach a sharply peaked Gaussian. In section 7.6 we estimate 
the dust spectrum width upper limit for which the dynamics we model dominate. In the 
absence of a detailed sticking and destruction model for collisions the extra complication 
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involved with using a spectrum does not guarantee extra precision for the present calculation. 
To encapsulate the unknown physics of dust-dust interactions, we presume that a collision 
between two dust grains has a probability p to result in sticking. If such a collision does not 
result in sticking the grains are presumed to be unchanged. Accordingly we can see that, if 
N is the dust-dust collision rate then ^f- = pNrrid and so using the equation for m d (r d ), 

drd Pat m 

As dust grain size is a key parameter whose change as a function of the distance from the 
star we investigate, we will require that dust grains do not significantly migrate radially. 



2.1. Equation of grain growth 

A simple "particle-in-a-box" collision model implies that 

N = — a dd v c = — a dd v c = — , (2) 

m d H d m d p d r d H d 

where H d is the dust disc scale height, p and S rf are the dust disc volume and surface 
densities, f2 is the orbital angular velocity and v c is the dust-dust collisional velocity. The 
dust and gas scale heights need not be equal, the latter being given by H g = c s f2 _1 , where c s 
is the gas thermal speed. We can also define a velocity v h = H d VL and write Eq. 2 in terms of 
a non-dimensional ratio v c /vh- Our method of handling turbulence lends itself more readily 
to Eq. 2 however. 

If the velocities of the dust grains are uncorrelated, then their collisional velocity v c is 
also their random velocity dispersion and so, akin to the formula for the gas scale height 
H g = c s /n, we have that dust scale height H d = v c /Q. If the ratio v c /H d in Eq. 2 is 
independent of grain size r d , then N ~ rj 1 . For v c = VtH d and p d = lgcm -3 the growth rate 
from Eq. 1, ^ ~ 107 (y^u) 3 ^ is independent of dust size and the formation of meter 
sized bodies will occur in ~ 1 year. As without turbulence, gravitational instability of the 
disc will occur for sub-centimeter sized grains (Eq. 74); this process occurs well within the 
time scales required by existing constraints (e.g. D'Alessio et al. (2005)) even considering 
settling times. 

If, however we allow the velocities of the dust grains to be correlated, as would occur 
for grains in the same eddy, there is no reason to presume that v c = VtH d or that v c and H d 
are even closely related. If the dust couples strongly to the gas and the turbulence causes 
significant mass mixing throught the gas disc, then the gas and dust will share a common 
scale height H d = H g . In addition, very small grains will have the same bulk velocity as any 
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turbulent eddy in which they reside, and will interact with collisional velocities corresponding 
to the dust thermal speed, i.e. v c = v t h- The gas sound speed c s is much greater than the 
thermal speed of even small dust grains v t h so the dust growth rate would be prohibatively 
low if such conditions held for large ranges of r d . 

We can, however, imagine that as the dust grains grow and couple less strongly with 
the gas, the turbulence can collide dust grains together with collisional speeds v c ~ v tur buient 
that are much greater than v t h- Eventually the growing dust grains will settle out of the 
gas disc, so that H d will decrease. For very large dust grains, with mean free paths (with 
respect to each other) longer than the largest scale eddies, the effects of the eddies will be 
uncorrelated, causing v c to equal QH d , drastically increasing the growth rate (Eq. 2) above 
the value for very small grains and speeding the onset of gravitational instability. 

To facilitate deriving the formulae for v c and H d in different regimes, we define = p d r d , 
which is proportional to a single dust grain's surface density. Then, presuming a constant 
Pd, we can combine (1) and (2) to obtain: 

the form of the growth equation we will use for the rest of the paper. The variable <f> will 
later help clarify the effects of different dust grain densities on time scales and critical grain 
sizes. As is the dust density in the disc, T, d jf^dt is the surface density of the dust sheet 
that a dust grain travels through in time dt. Hence (3) states that the sticking-parameter 
modified surface density of that sheet can be added to the penetrating dust grain's surface 
density, much like a ball moving through cling wrap. 

The basic question we seek to answer is whether the grain growth is sufficiently rapid 
that the gravitationally unstable stage of planet formation can be reached long before the 
total time scale available for planet formation (~ 10 6 yr, D'Alessio et al. (2005)). In our 
study, we need to consider that very rapid collisional velocities could be too high for sticking 
to occur and instead result in the destruction of the grains. We also need to determine 
how large the grains must become for them to settle out and allow gravitational instabilities 
to occur. While there exists a dust scale above which the effect of turbulence on the dust 
velocity dispersion becomes sufficiently weak for the dust disc to be gravitationally unstable, 
that size scale will be far greater than the initial size of the dust, as well as that of the critical 
size in the absence of turbulence. A last consideration is our requirement that the growth 
occurs before significant orbital migration. 

To address these issues, we determine various regimes for v c and H d before the onset 
of gravitational instability. These regimes depend on the quantities r d , c s , the turbulent 
viscosity parameter a and the mean free path of the gas l m f P - We then use (3) and the 
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values for the ratio of v c and H d to calculate the peak collisional velocities and the total 
growth time before the onset of gravitational instability. 



2.2. Disc model 

We assume the disc has an essentially isotropic turbulent viscosity 

v T ~ %A M ~ ac s H g , (4) 

where % and Am are the maximum turbulent velocity and length scales in our assumed Kol- 
mogorov spectrum and a is the Shakura-Sunyaev viscosity parameter (Shakura & Sunyaev 
1973). Using Xm/vm ~ 1/^, H g VL = c s , and (4) it follows that 

A M ^ a*H g (5) 

v M ^ a^c s . (6) 

If a is extremely low, the GW process will proceed as in the absence of turbulence. It is likely 
that at least dust-gas shear (e.g. Ishitsu & Sekiya (2003)) prevents a from dropping to values 
where it could otherwise be ignored. Our model will involve a progression of dust growth 
regimes whose order requires approximately a > 2 x 10~ 6 . Our model needs modification for 
values of a below that lower bound and if the turbulence is adequately weaker yet that the 
dust settling time becomes significant then our approach fails completely. Thus although we 
are conceptually extending the a = model of GW, our present techniques do not apply in 
the turbulence-free limit. 

We take a minimum mass solar nebula (MMSN) model with the following scalings from 
Sano et al. (2000): 

Z 9 (R) = S 9 o (^j) 2 (7) 

= do) 

where S s0 ~ 1.7 x 10 3 gcm' 2 and T = 280K. We presume that S d = S s x 10" 2 , and 
cjgg = 10~ 15 cm 2 where a gg is the neutral gas-gas cross-section. The average gas molecule 
mass m g = ixtuh where \i = 2.34 is the mean molecular weight of the gas. The above values 
are presumed to remain constant in time. 
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2.3. Incorporating the role of disc turbulence and drag 

As the dust grains grow, the effect of turbulence on their collisional velocity and scale 
height evolves. For dust grains smaller than the mean free path of the gas, that is r d < 
l m fp = ~~ ; the gas-dust interaction is characterized by Epstein drag, whereas for r d > l m f p 
we have Stokes drag (Landau & Lifshitz 1959). Throughout this paper we will use E and 
S as subscripts refering to behaviour with Epstein and Stokes drag respectively. For a dust 
grain interacting with a turbulent eddy of size A and relative speed v = v(0) at some initial 
time t — 0, the time evolution of v subject to Epstein drag is given by 

dv = _ 2 P E ^c ±v = _^ (n) 
dt m d t e 



while for Stokes drag 



dv _ 67ir d c s u ^ _ 2irl m f p p g r d c s ^ _ v ^ 
dt m d m d r s ' 



where p g is the gas density, a gd = nr\ is the gas-dust cross-section and v = pgl ™ fp chosen so 
that the drag equations are equal at the boundary between the Epstein and Stokes regimes 
r d = l m fp- We can combine the above two equations by writing 



where r = Mqx(te, ts) and 



and 

We note that even though the gas is turbulent, the scale of the smallest turbulent eddy is 
generally larger than the scale of the dust grains during the growth phases we consider. This 
is why we do not use the high Reynolds number form (Landau & Lifshitz (1959)) of the 
Stokes drag (which varies as v 2 instead of as v). This will be justified a posteriori later. 

Eqs. (11) and (12) have the same form whose solution is 

v(t) = v(0)e~r . (16) 



where again r = Max(r s , r E ) from (14) or (15). If v(0) = V\, the turbulent velocity for an 
eddy of scale A, Eq. (16) can be integrated to show that the dust grain will not exit the eddy 
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in a time less than the eddy turnover, or destruction, time t\. Upon being released from the 
eddy at its destruction, the grain's speed will have changed by 



For a given stopping time r, Av from the above equation has its maximum for the largest 
scale eddies and so 



is the dust's random velocity with respect to the Keplerian disc and the latter equality follows 
from Eqs 5 and 6 

3. Calculating the collisonal velocities and scale heights v c and H d 

To facilitate incorporating the calculations of the previous sections into this and the 
following sections, we have collected key variables used in this paper into Table 1. 

VMRJ developed an elegant model, later improved by MMV, describing the velocities 
of dust grains in disks with Komolgorov turbulent spectra. Dubrulle et al. (1995) used their 
results to calculate dust scale heights in such discs. Unfortunately, VMRJ made recourse 
to an auto-correlation function (ACF) that was shown by MMV to behave poorly for small 
time differences (where the most interesting physics occurs). Neither MMV (their equation 
5) nor Cuzzi & Hogan (2003) strongly motivate their different choices of ACF and we show 
through direct physical argument that some of their results (and those of Dubrulle et al. 
(1995)) are unphysical, raising doubt about the validity of their ACF. The question of what 
ACF reproduces our results, if any, is beyond the scope of the present paper. 

The need for an ACF stems from VMRJ's consideration of the question of a spectrum of 
dust grains interacting with turbulence. We instead make the simplifying approximation that 
all the grains are the same size r&. For our calculations we note that grains exiting an eddy 
will promptly enter another eddy of comparable size. This assumption is straightforwardly 
justified if the turbulence is characterized by smaller scale eddies embedded fractally within 
larger scale eddies. If instead smaller scale eddies surround larger eddies, then justification 
of the above assumption follows from the Kolmogorov hypothesis that eddies of a given 
size cascade by interaction with eddies of comparable size and the fact that smaller eddies 
have lower speeds and less influence on a dust grain speed than a larger eddy. As shown 
schematically in Fig. 1, there will be large eddies, whose effects on the dust grains are 
too correlated to contribute significantly to the collisional velocity, small eddies too weak to 




(17) 




(18) 
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significantly effect the dust grains and a "Goldilocks" scale of eddies "just right" to dominate 
the interaction of dust grains. By making the simplifying assumption that all the grains have 
the same size, it is a simple task to identify the Goldilocks eddy scale as a function of grain 
size and calculate its effects. 



3.1. Collisional Velocity of Grains 

To estimate the characteristic collisional velocity between two dust grains we assume 
that the grains' individual speeds are each Av from (17), determined by having exited 
separate eddies of comparable size. We then estimate the collisional velocity v c as the 
average relative speed of dust grains from head-on and catch-up collisions: 

u c (0) ~ ^(Au(0) + Av(0)) + ^(Au(0) - Au(0)) = Au(0) =v x (l- e"^) , (19) 

where the latter equality follows from (17). Since Av represents a linear combination of two 
dust grain speeds, and both sides of (11) and (12) above are linear, we then know from (16) 
and (19) that v c evolves as 



v c (t) = v c (0)e—r = v x [1 - er J e" . (20) 
Averaging (20) over an eddy turnover time from t — to t\ then gives 

vc(x) = ±- v e (t)dt = ^ (1 - e-*) 2 , (21) 
t\ Jo x 



where x = ^. For a Komolgorov turbulent spectrum, v\ = % j 3 and A = V\t\ which 
we can use with (5) and (6) to write: 

v x = {tVLx)^ sfac s . (22) 

Together, Eqs 21 and 22 give 

i fl —x\^ 

v c (x) = (rfi)2 1 ~ 6 ' V^c s . (23) 
v x 

The Goldilocks eddy scale is the one that maximises v c (x). Equations 24, 25, 26 and 27 below 
describe the dust-dust collisional velocity v c in four different size regimes and are gathered 
in Table 2. 
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3.1.1. Maximuma of v c (x) 

From Eq. 23, maximizing v c (x) gives 

Max[v c (x)] = v C:2 = BVrVL^fac s (24) 

with B = 0.53 and the maximum occurs at x — x max = 2.3. As long as there exists a 
A for which x = x max we will use Eq. 24 as the turbulence induced collisional velocity. 
However, a Komolgorov turbulent spectrum has both a maximum and a minimum scale 
such that t M > t x > t m where t M = = ^ is the turnover time for the largest scale eddy 

and t m = \J 1™* p q is that for the smallest. We seperately treat those dust sizes for which 
Goldilocks eddy scales do not exist. 



3.1.2. Small r cases 



For small enough r, the minimal value of x = ^ > x max and we cannot use Eq. 24. 
However for such large x, (1 — e~ x ) 2 is close to 1 and we can approximate (23) as: 



i . i 



v c (x) = v cA {x) = [ — ) V®c s = f — J Ty/ac s . (25) 

This approximation is useful because the time scale of the regime (regime 1.1) where we will 
use it is only weakly dependent on its final size as can be seen in Table 4. 

For sufficiently small values of r, the values of v c from (25) will be less even than the 
dust's thermal speed. For such grains we then just use the dust thermal velocity for v c , 



<V = r,„(o) = ,v,,(o) = r.^/^ = A /^f . (20) 




where m g is the molecular weight of the gas. 



3.1.3. Large r case 

For large enough r, the maximal value of x = tj ^- = ^ < x max and again we cannot use 
(24). In this regime we have 

v c = v c , 3 = v M - (1 - e~ x ) = c s ^faVLT \1 - e~^J , (27) 
as that maximizes Eq. 23 given the constraints on x. 
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3.1. 4. Maximal v c (t) 

Equations (25), (24) and (27) allow us to maximise v c with respect to r. We find 

Max[v c (r)] = v c , T = C^fac s , (28) 

for C = 0.407 at x — 1.26, which is less than x max . Accordingly, the maximum v c (t) occurs 
in the regime for which v c = v c ^. 



3.2. Scale height H d 

In addition to affecting the collisional velocities, turbulence also affects the dust scale 
height Hd by preventing settling. We can use Eq. 13 to calculate the settling time of a 
dust grain in the disc. Ignoring turbulence for the moment, a grain at height z above the 
midplane moving towards the midplane with speed v = % feels a vertical force 



dt 

d 2 z GM z Idz 



dt 2 R 2 R rdt 
which gives a terminal velocity (^# = 0) of 

dz GM 



(29) 



TZ = -n 2 TZ. (30) 



dt R 3 

Accordingly, z(t) = z(0)e~~^ , where the settling time t s = (fi 2 r) \ 

In the presence of turbulence however, there will be "updrafts" of turbulent eddies 
supporting the grains against the gravity. Balancing the force of gravity with the force of a 
largest scale eddy on a grain at rest gives: 

Q 2 z = V JL = (31) 

T T 

which we rewrite as 

/5 



- Or*' (32) 
As long as ->/a/(fh~) < 1, force-balance will occur for z < H g and hence we can use z from 
(32) as Hd- For smaller grains however, implicit assumptions of isotropy of turbulence break 
down at the boundaries of the gas disc and so we limit the scale height of the dust to that of 
the gas Hd = H g . Accordingly we have that for y/a/iVLr) > 1, the dust does not settle and: 



Hd,i = H g , 



(33) 
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while for y/a/(Qr) < 1 the dust settles and instead: 




The results are gathered in Table 2. It is worth noting that the above is well approximated by 
the distance the largest scale eddies (which have the greatest effect on the dust's movement 
with respect to the Keplerian disc) can carry a grain in a single settling time, bounded above 
by the scale height of the gas disc and below by the distance a single largest scale eddy can 
carry a grain. 

In our calculations of H d 2 we have implicitly assumed that the grain growth time scale 
is longer than the settling time scale. If we calculate the growth time scale | and the settling 

time t s at the onset of settling (r = ^), for Epstein drag we find that the settling time is 

_ i 

shorter by a factor of 2pB ° 1 |K For a > 1CT 8 then, our assumption is justified. 



3.3. Comparison to previous work 

MMV improved on the work of VMRJ by noting that the cut-off of turbulence at small 
scales effects the dust's velocities, and found all the regimes we have (though they ignored 
the possibility of thermal interactions). However, the behaviour they found in the regime 
that corresponds to our section 3.1.2 differs from our results. In this regime, the velocities 
of dust grains inside a given eddy will be extremely correlated because the grains are very 
small. It follows that their interactions will be purely thermal except inside the thin "skin" 
layer of each eddy into which dust grains can penetrate before losing any initial velocity 
relative to the eddy. We can approximate this by sending a test grain into such an eddy 
with initial velocity v. The test grain will interact with those dust grains already trapped by 
that eddy with its initial velocity for its stopping time r before decelerating and becoming 
trapped. As the grain will remain trapped for an eddy turnover time t\, the grain's average 
collisional velocity will be: 

vr + x (t x - r) vr 

iJaverage ~ , — , • \°°) 

The above physical estimate not only scales with r as (25) but the two equations are in fact 
equal. Because the results of MMV show a different and unphysical scaling in this regime, 
we find it likely that their ACF and its implications are unphysical. 

As the approach of Kornet et al. (2001) and Stepinski & Valageas (1997) (using the 
results of Cuzzi et al. (1993)) involves averaging over the turbulent energy spectrum, their 
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results are only weakly dependent on the small scale turbulence and they do not consider 
this extremely correlated regime just discussed for the collisional velocity. It is worth noting 
that Cuzzi et al. (1993) uses a non-Kolmogorov turbulent energy spectrum (Meek & Jones 
1973). 

Another difference between our work and that of VMRJ or Kornet et al. (2001) is in the 
behaviour of large (tQ > 1) dust grains, for which settling occurs. Our calculation of the 
dust scale height differs from Dubrulle et al. (1995) (using the work of VMRJ) and Kornet et 
al. (2001) who find the dust scale height behaving as while ours behaves as r -1 . As we 
derived our scale height from direct considerations of force-balance however, we feel justified 
in using our results. 

4. Regimes of the equation for grain growth 

In the previous sections we have derived a progression of four formulae for v c and 
two for Hd that apply as dust grains grow, as gathered in Table 2. As long as the dust 
grows monotonically, for any values of R and a we have a progression of up to five regimes 
for the ratio and hence Eq. 3. For R and a in the allowed range of this paper ( 
2 x 10~ 6 < a < 10~ 2 and 0.5 AU < R < 8AU), Table 2 summarizes the regimes which we 
label by i.j where % corresponds to v c ^ and j corresponds to H^j. Table 2 can be used to 
calculate the grain sizes for which regime transition occurs, and to solve Eq. 3 for the time 
it takes for the grains to grow through each regime into the next. As these calculations are 
straightforward, tedious, and repeated for 5 regimes, they are collected in the Appendix 1. 
Tables 3 and 4 contain the results of these calculations. 

As our time scale calculations involve ratios of v c and Hd we can see that once the grains 
begin to settle (regimes 2.2 and 3.2), the factors of a will cancel in Eq. (3), resulting in the 
lack of a dependence of the time scales seen in Table 4. However, as seen in Table 3, the 
grain size range over which the regimes occur still depends on a because the beginning and 
conditions for each regime depend on a. Thus the time scale for each regime, and the overall 
growth time to the gravitational instability phase still depends on a. 
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5. Results for growth time scales and sizes 



5.1. Grain sizes at the end of each regime 



Figure 2 shows the grain sizes calculated from the formulae for 0j for the end of each 
regime of Table 3 as a function of disc location. The calculations use a fixed p d = lgcm" 3 , 
so r d = r d (0) oc . The regimes are labeled in the Figure captions. The solid line is the 
curve pdlmfp, which represents the boundary between Stokes (left) and Epstein (right) drag 
regimes. 

An important feature of Fig. 2 is that there is a maximum grain size in each regime 
which falls along the line bounding the Stokes and Epstein drag regimes. This arises because 
all of the equations for <pf of the previous sections that define the boundaries between regimes 
are increasing functions of R for Stokes drag and decreasing functions of R for Epstein drag. 
Because the drag viscosities are equal at the boundary between Stokes and Epstein regimes 
(Sec. 2.3), for each regime of Sec. 4, there is a triple point in the plots of 0/ vs R where 
4>e = <Ps = Pdlmfp- The triple point is formed by the intersection of the curves corresponding 
to Epstein and Stokes drags, and the line which represents this boundary as a function of R. 
Compared to the Stokes regime, the Epstein regime curve is less dependent on the various 
parameters such as pd and l m f p . 

The triple point of the curve r = Q^ 1 is particularly significant because that curve 
defines the grain size scale for which grain infall from the drag on the gas is important 
(discussed in Sec. 7.2). That this is where the infall is maximized is not surprising as it 
occurs when all of our time scales, the stopping time r, the orbital time Q^ 1 and the settling 
time t s are equal. For our disc model this triple point is independent of a and is given by 



Although we first present results without including dust migration, note that due to 
the monotonic behaviour of the boundaries in a given drag regime, inward migration would 
cause grains of a constant size to rise above the size curves of Figs. 1 and 2 in the Stokes 
regime but fall below the curves in the Epstein regime. As a grain grows from the smallest 
to the largest within a regime, it will experience Epstein drag for R above the triple point 
of the regime, Stokes drag if R is less than the triple point of the prior regime, and it will 
transit from Epstein to Stokes drag if R lies between the two triple points. Figure 2 shows 
that because of this behaviour, the peak grain size occurs for 5 < R < 8AU for dust grains of 
r d ~ several meters. This is the size at which gravitational instability takes over and marks 
the transition to growth regimes subsequent to those studied in this paper. 
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4 




(36) 
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5.2. Total time scales to grow to gravitationally unstable regime 

In figures 3, 4 and 5 we plot the total time scale tr to grow through all regimes 0.1 
through 3.2 and thus the total time to reach the gravitationally unstable phase as calculated 
from Tables 3 and 4. A key feature to note is the a dependence of tr, as exemplified in Fig. 
3. For regimes 0.1 through 2.2, the time spent in a given regime is a decreasing function of 
a. The duration of regime 3.2 is an increasing function of a. If follows that for any given 
R there will be an a that minimizes tr- For R = 3AU this value is a ~ 3 x 10~ 2 while 
for R = 8AU it is a ~ 10~ 3 . Provided a and R are within the range described in Sec. 
7.1, this allows us to determine the minimum and maximum values of a for which the total 
growth through all regimes is less than the observationally constrained time scale (< 10 6 yr) 
of planetesimal formation. 

Figure 5 is analagous to Fig. 3 but for different dust grain material densities and 
dust /gas surface densities in the disc (see Fig. caption). 

Except for regime 0.1, equation 14 and the Epstein regime time scale equations (53, 
58 and similar) depend neither on pd nor l m f p . As our regime boundary and maximum 
velocity equations are cast in terms of r, the use of <fi therefore makes clear the Epstein 
regime's lack of dependance on p d and l m f p . However, the mass and radius of the grains 
vary as Pd~ 2 <p 3 or p^ 1 at various stages of the evolution so the actual grain radii at which the 
regime transitions occur do vary. In the Stokes regime, density plays a non-trivial role in the 
effects of turbulence on the grain trajectory (15) and therefore in determining the boundary 
between the two drag regimes as well as the details of the regimes of jf^. In figure 5 we plot 
the time scales for p d = 0.5gcm~ 3 . 

From (3), we see that increasing the dust disc surface density will decrease the growth 
time scale. However, in our equations frequently appears in a ratio with T, g . The effect of 
changing E d while keeping S s constant is greater than that of changing both equally, as can 
be seen in Fig. (5). The total surface density primarily affects the condition for gravitational 
instability which determines the time spent in regime 3.3. This dominates the growth time 
scale at large R. Changing the dust surface density alone has a strong influence throughout 
the disc. 



5.3. Dead zones reduce the total growth time 

The MRI (Balbus & Hawley 1991) instability may be the primary source of turbulence in 
a protoplanetary disc but the entire gas disc need not be sufficiently ionized for the instability 
to occur. The MRI-stable region of the disc is labeled a dead zone (Matsumura & Pudritz 
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(2005), Sano et al. (2000)) and will have a smaller value for a than the live zone. If there 
is a dead zone for approximately the gas scale height over some range of R, then we can 
simply treat the live and dead zones as having seperate values of a and merging them in 
our radially depedendent calculations, while noting that dust unable to decouple from the 
turbulence on a viscous time scale could be deposited at the outer edge of the dead zone. 

If instead we consider a turbulent disc with a dead zone near the midplane of the 
disc (with a height Hp) for some range of R (but live for a significant thickness above the 
midplane) we can approximate the growth time scales by calculating the time scales for the 
live zone turbulence until the dust disc is contained within the dead zone Ha = Ho at which 
point we proceed with the dead zone turbulence. Noting that the time scales of regimes 2.2 
and 3.2 depend on a only for their beginning and end conditions, respectively, and that the 
onset of regime 2.2 is precisely the same as that of settling it follows that the dead zone a 
only effects the time scale of regime 3.2 and so the effect is simply calculated. Note that 
halfway settling the dust, H d = 0.5H g occurs for roughly 

within regime 2.2 (Table 2). The time scale behaviour with a dead zone is important because 
the duration of regime 3.2 is an increasing function of a, whereas regimes 0.1 through 2.2 
are decreasing functions of a. Regime 2.2's time scale, as discussed above, will only depend 
on the live zone a. Accordingly, as the contrast between the values of a in the live and dead 
zones increase, the the total time for planetesimal growth decreases. In particular therefore, 
having a dead zone, as opposed to no dead zone, decreases the total time below that of a 
disc with only one of the two values of a. 

The role of a dead zone in lowering the growth time scale is shown in Fig. (4). 



6. Constraints on the range of a, p and R for which gravitational instability 

occurs 

Armed with equations for growth rates and collisional velocities of the dust grains, we 
can identify the range of values of a and p that will permit planetesimal growth in a realistic 
turbulent disc on observable time scales for ranges of R. 

While we are restricted by our regime definitions and regime ordering to values of 
a > 10~ 6 , the effect of a on the estimated time scales for planetesimal growth is pronounced. 
We now discuss the effect of a and p on the range of values that allow for planetesimal 
formation on the observationally relevant time scales. 
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6.1. Constraints implied by the maximum v c 

We can place a maximum value for a that still allows planetesimal formation by studying 
its influence on the the maximal collisional velocity. We have used only the simple sticking 
parameter c to characterize the extent of sticking between interacting dust grains. However, 
for large enough values of v c , the dust grains will destroy each other on impact. A more 
complete dust grain interaction model would allow us to constrain v c and thus a because the 
maximum v c varies as ^Ja. from (28). For a disc with a dead zone, the maximum collisional 
velocity in the dead zone will be lower than that in the live zone because of the lower value 
of a. This helps facilitate planet formation even if the live zone a is otherwise too large. We 
quantify this below. 

If the dust scale height for r = fi -1 approximately defines the onset of the dead zone 
(Ha = yfaHg « H g for either a) then we can find the maximum live zone collisional 
velocity from (28): 

v c ,max = V®c s (l - e^ 1 ) ~ 0Ay/ac 8 . (38) 

If instead the end of regime 2.2 defines the beginning of the dead zone, then from Eqs 
24 and 65 we have i 

Bc s ~ 0.35x/ac s , (39) 

and a dead zone height from Eqs 34 and 65 of 

H D = ^fax max H g ~ 2^/aHg. (40) 

If the height of the dead zone is roughly half that of the gas disc, then Eq. 37 gives the 
condition for the dead zone height to equal the dust scale height. The maximum live zone 
collisional velocity occurs for that value of r and is 

v c ,max = V2a*Bc s ~ 0.75«t Cs . (41) 



We plot results for live zone values of a = 10~ 3 and a = 10~ 4 in Fig. (6). The 
maximum collisional velocities occur for dust grains of centimeter to meter scales. The top 
solid line in each panel represents the maximum collisional velocities without a dead zone 
and the bottom solid line represents the maximal live zone collisional velocities for a dead 
zone with Hp = 0.5H g . From common experience, we infer that it is unlikely for dust grains 
colliding with speeds greater than 30km/h to stick from the collision (and more likely that 
they obliterate). This constraint is consistent with Youdin (2004) who finds limits for v c 
near 5ms _1 ~ 18km/h from the physics of dust grain collisions. Therefore we interpret Fig. 
(6) to imply that, while for a live zone a near or above 10~ 3 , a significant dead zone (with 
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a lower a ~ 10~ 4 ) is needed in order to avoid the collisional destruction and enable planet 
formation in the absence of additional processes like vortices. However, if the live zone a 
were as low as 1CT 4 then a dead zone would be unecessary. 

6.2. Constraints implied by the total growth time scale through regime 3.2 

The discovery of a planet younger than a million years old (D'Alessio et al. 2005) 
requires, at minimum, that planetesimals grow to the gravitational unstable phase on a time 
scale substantially shorter. For R < 3AU and p — 1 our total growth time scale to reach this 
phase, tj 1 , satisfies tx < 10 4 yr for 10~ 3 < a < 10~ 4 , allowing at least an order of magnitude 
leeway while still staying well below 10 6 yr (Fig. 3). We also see that as a falls below 10~ 4 
the total time scale rapidly rises. This allows us to place a lower bound on a greater than 
the bound required by the validity (section 4.2) of our regime ordering in Table (2). 

7. On the validity of key approximations 

7.1. Range of a for which Table 2 is applicable. 

Both the existence and order of the regimes for v c and Hd (Table 2) can vary based on 
the strength of the turbulence (a). For example, if a — then the only form of v c of ours 
that applies is thermal (Eq. 26), and our derivation of H d does not apply. 

For 0.5AU < R < 8AU no changes from the canonical regime ordering occur for 2 x 
10~ 6 < a < 1/3. As a drops below about 10" 6 our present regime ordering of Table 2 breaks 
down as shown below, though a new ordering (and possibly new regimes) could be found 
through a similar approach, as long as the settling time remains short (section 3.2). 

If a < x^ ax , then ^ < — ^ — . It follows that the transition from Vh,i to Vh,2 will occur 
before that from f Cj2 to v c ^. For large grains, a possible deviation from the bottom-to-top 
canonical ordering shown in Table 2 arises if gravitational instability is triggered before i> C;3 , 
which for 0.5AU < R < 8AU requires a < 10 -8 . For small grains, possible deviations from 
the canonical regime ordering occurs when = ^ (Sec. 4.3) and when equations 51, 52 
(see appendix) hold for r > t m x max . In this case, the weakness of the turbulence delays the 
transition from thermal collisions to turbulent collisions. 
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7.2. Constraints from comparing headwind induced infall and grain growth 

time scales 

We have made the approximation that and T, g (R) are time independent and that 

dust grains do not migrate. However, while gas depleted by accretion is replenished from 
the outer regions of the disc or envelope, dust migration cannot be so simply ignored. The 
gas in a protoplanetary disc is partially pressure supported and so rotates at a sub-keplerian 
velocity. The dust grains will then feel a headwind and begin to spiral in towards the star 
(Weidenschilling (1977), Nakagawa et al. (1986)). As seen in Nakagawa et al. (1986), this 
sub-keplerian rotation behaves as their eq 1.9 (which omits a fraction sign): 

v rb,g = (1 -r])v K , (42) 
where v k is the Keplerian velocity and 

R 2 dP _ 13 / c_s \ 2 , . 

V 2GMp g dR 8 \QRJ { } 

parameterises the difference between the gas orbit speed and the Keplerian speed. The gas 
disc also accretes viscously, with time scale t v = a ^ 2 H and accretion velocity v v ~ If we 
assume p g » p, our r _1 corresponds to the D of Nakagawa et al. (1986), so adding in the 
viscous radial velocity we rewrite their eq 2.11 for the infall speed as 

T OL(? I 

K -- 2 " tfR TWsrr^-" +t » (44) 

where v r = 2nQ 2 R-— ! — TI is due to the headwind and v v = * , 2 is due to the viscous 

' l+(S2r) till l+(S2r) 

accretion flow. As the infall due to viscous accretion flow v v < 1.2 x 10 -5 (j^) -jf is too 
small to allow for significant infall within 10 5 years we neglect the term. The maximum 
Vr(t) then is r]VtR ~ |^ and occurs, for any R, at r = f] -1 . For small r, Vr scales as r 
while for large r, Vr scales as r _1 . 

Migration can be ignored only if the decceleration time scale of the infall speed —Vr/Vr 
is less than the infall time scale R/Vr for the most affected grains. For r > Q^ 1 , Vr oc t and 
the decceleration time scale can be written as r/f for any drag regime. We can then cast 
the condition to ignore migration as the requirement that the ratio Rjy R < 1. In Fig. 7 we 
plot 1/c times this ratio. The jump in the plot at R ~ 6AU is due to the fact that ^ = ^ 

but ^ — i and R = 6AU is the triple point of r = Q^ 1 (sec. 5.1) where the drag regime for 
the appropriately sized grains changes (as seen in Fig. 2). In Fig. 8 we plot the R inside of 
which the headwind infall can be ignored, as a function of p. It follows that infalling grains 
will pile up at this radius, leading to a disc dust density enhancement and opportunity for 
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enhanced planetesimal formation. This in turn, allows more flexibility in the value range for 
p that still facilitates rapid planetesimal formation in the absence of other physics. 

There is a subtlety however. As we can see in Fig. (2), inside of R tp (sec. 5.1, for our 
disc R = 6AU), the curve r = (second curve from the top in the right panel) lies within 
the Stokes drag regime, while for R > R tp Epstein drag applies. In the Stokes regime, the 
curve is an increasing function of R so as grains are transported inward, their r will increase 
above and grain growth past the curve will always occur. For the Epstein drag regime 
however, the inward transport will keep the grains below the r = Q^ 1 curve; here if grain 
growth at a given radius is slower than the inward transport, the grain will never transit the 
curve, preventing further growth and forbidding the onset of gravitational instability. 

To see whether this circumstance arises, we calculate the rate of change of <p c (where 
<p c = -if- is the critical size required to satisfy te = f2 ) due to the radial transport and 
compare it to ^ from Eq. (3). We find 

dcj> c rf3S s 3<j) c dR 3 117 £ g Qff g 2 

It = dt-T = ~2R^ = 2^ QR = -36^- (45) 

— - 3 

where £ 9 scales as R~? from (7) and 

^ = pZ d Q (l - e~ l ) (46) 

for grains with = <p c . For the disc model we have chosen in Sec. 2.2., we find for = C 
that 

<**/* Kpl . g f * yV5#y\ (47) 



d<p c /dt * ' VAuy v 100 

In that portion of the disk where this quantity is less than 1 in the Epstein drag regime, 
the inward transport prevents the grains from growing above r = Q' 1 . For our canonical 
disc parameters this occurs everywhere that r = Q^ 1 occurs for Epstein drag (i.e. R > Rt p ). 
Because the time scale to reach C is an increasing function of R while the total mass 
contained within an annulus of width dR is a decreasing function of R, changes in as 
the dust falls inward do not alter the basic conclusion. Increasing S 9 increases the triple 
point Rt p of (Eq. (36)) and increasing the ratio S^/S. decreases the radius where (47) drops 
below unity. Thus the critical radius above which planetesimal growth is prevented is indeed 
R tp . The effect remains important if the basic disc parameters are altered, but the radius of 
the triple point and and the critical radius where (47) drops below unity and the triple point 
radius (36) can change. 

The important upshot of the effect just described is that gravitational instability cannot 
be easily reached at radii larger than the triple point (discussed in section 5.1) of the r = 
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curve without additional physics like vortices, or a significantly different disc model. It 
is interesting, if not unsurprising in hindsight, that the drag induced grain agglomeration 
benefit from vortices (Barge & Sommeria 1995) is maximized for the same r = Q^ 1 condition 
that gives the maximal infall speed we found from (44) . 

The conclusions drawn from (47) and the limits on p from Fig. (8) are essentially 
independent of a because its effects on v c and in Eq. (3) for the relevant regimes cancel. 
This allows us to use Fig. (8) to constrain our sticking parameter, indepedent of a (at least 
for R < 8AU). This makes our constrains on a more meaningful. 

7.3. Neglect of gravitational enhancement to cross section 

Our model ends at the point when gravitational instabilities take over. We have calcu- 
lated whether gravitational enhancement to dust-dust collisions plays any role before that 
time. We find the effect to be small (typically under 5%) for all but the largest a even for 
the largest scale dust. 

7.4. Use of linear vs. quadratic in v drag forces 

Even though the present flow is turbulent, the laminar approximation for the Stokes 
drag was used because the grains were presumed to be smaller than the smallest turbulent 
eddies at all times. Therefore, unlike a macroscopic object in a turbulent flow, the grain 
would see an effectively laminar flow. 

This approximation can be justified a posteriori: It turns out that only very close to the 
star, do the largest scale dust grains just before gravitational instablity become comparable 
to the smallest scale turbulence. We neglect the change in drag law that this would imply 
because those eddies both are far less massive than the grains and are traveling at velocities 
nearly equal to these dust grains at the onset gravitational instability and so will have 
negligeable effects on the grain trajectories. 

7.5. Sticking 

The most uncertain aspect of our model lies in the sticking mechanism. A value of p = 1, 
which represences 100% efficiency for collisions resulting in sticking, seems unlikely. Consid- 
ering p to be a gestalt parameter combining both all the possible sticking and destruction 
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results of collisions and examining Eq. 3 we can see that our time scale T depends inversely 
on p. Comparison with sec. 7.2 we see that decreasing c strongly decreases the outer value 
for R allowed by the headwind induced infall, although this will only become important if 
p nears 10~ 2 . If we had a well developed sticking model we could calculate the maximum 
collisional velocity before grain collisions result in destruction. Accordingly, our ability to 
constrain a, either through disallowing grain destruction through high velocity collisions or 
through time scale arguments is limited by the quality of the sticking model. 



7.6. Dust size spectrum 

In our model we assume that all the dust grains share a common radius r^, and in 
section 2 we justify this assumption by appealing to an eventual Gaussian spectrum. As 
we saw in section 7.2, the differential rotation of the dust and gas discs caused by partial 
pressure support of the latter induces size dependent velocities in the former, the radial one 
given by Eq. 44 and an orbital one (Nakagawa et al. (1986)): 

-rfw- (48) 

An adequately broad dust size spectrum would then affect the grain-grain interactions by 
changing the form of v c . We can evaluate the importance of this effect by calculating the 
width of the dust spectrum (Ar/r) required for the difference in headwind induced velocity 
to equal the velocity induced by the eddies which dominate dust-dust collisions in our model 
from Eq. 17. We therefore calculate the minimum Ar/r for which a spectrum would 
influence our results. The larger the Ar/r we calculate, the broader the spectrum our model 
can accommodate. For 0.5 AU < R < 8 AU and 10~ 4 < a < 10~ 2 , we find that the minimum 
Ar/r for which a spectrum would affect our results is 0.12, in regime 3.2. For Epstein drag, 
this is also Ar^/r^ while for Stokes drag this implies Ar^/r^ = 0.06. 

Draine (2006) gives an observationally derived spectrum for small grains that deviates 
from our Gaussian. However, for most of the planetesimal growth, the grains sizes that we 
model are much larger than the micron grain sizes in the Draine spectrum. While larger scale 
grain-grain collisions will produce some detritus down to the micron scales, we cease to model 
the spectrum on these scales once the bulk of the dust mass is predominately comprised of 
larger grains. There are few observational constraints on the grain spectrum on scales much 
larger than microns. This latter circumstance justifies our having pursued the simple "dust 
of a single size at any one time" approach herein. The benefit is that the basic physics can 
be identified at each stage of growth, the formalism is simple, and meaningful constraints in 
growth times can be calculated analytically. 
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8. Conclusion 

We have developed a simple model to estimate the effects of disc turbulence on plan- 
etesimal growth We show that turbulence has two competing influences: (1) It increases 
the collisional velocity of dust grains which increases their interaction for growth. (2) It 
increases the velocity which determines the scale height of the dust layer in the disc which 
decreases the volume density lowering the interaction for growth. Since the growth rate of 
grains depends on the ratio of these two velocities, the evolution of protoplanetesimal for- 
mation depends on the evolution of this ratio. Because the coupling of the gas to the dust 
evolves with grain size, the ratio is not a constant. Our work differs from previous studies 
by explicitly incorporting the dynamical evolution of dust grain size and the dust settling 
and considering discs with and without dead zones. 

We identify the range of turbulent strength, measured by a, the range of dust sticking 
efficiency, measured by c, and the disc radius range for which protoplanetesimal formation 
progresses to the gravitational instability phase on observtionally constrained time scales (< 
10 6 ) yr, D'Alessio et al. (2005)). The allowed ranges emerge from constraints on (1) the total 
time scale for planetesimal formation, (2) the infall time scale of the dust grains compared 
to their growth time and (3) the destruction of dust grains from high velocity collisions. 
We find that dust grain infall due to headwind gas drag prohibits reaching the gravitational 
instability phase outside a limiting radius that depends on the transition between Stokes 
and Epstein drag (the r = Q^ 1 triple point from sec. 5.1). For our disc model this occurs 
at R = R tp given by Eq. (36). It is possible that material piles up at R = R tp which may 
enhance planetesimal growth at this radius. This latter effect requires further work. 

For R < R tp , we find that for p > 10~ 2 and a < 10~ 3 turbulence does not catastroph- 
ically slow the needed planetesimal growth, as seen in Fig. (3), even though the time and 
size scales at which the growth occurs still depend strongly on those parameters. These 
results depend on the unknown physics of grain-grain interactions determining p, but the as- 
sumption that grain destruction results from collisions with velocities greater than 30km/hr 
strictly requires a < 10~ 2 for all discs if growth is to be fast enough. As determined in Sec. 
6.1.1, values of a ~ 10~ 3 for discs with substantial dead zones and a ~ 10~ 4 are acceptable. 

The condition that the grain infall time be longer than the grain growth time (as dis- 
cussed in section 7.2) also leads to an a independent minimum acceptable p as seen in Fig. 
(8). We find that as p falls below ^ the distance from the star at which gravitational in- 
stability can be reached and planetesimal formation occurs falls below 3 AU. For p = 10~ 2 
planetesimal formation won't occur outside of roughly 1 AU. Values of c below 10~ 2 are too 
low to produce rapid enough planetesimal growth in the absence of additional physics. 
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Although more detailed calculations with grain size spectra, collisional velocity spectra, 
and a more detailed grain-grain collision model are needed, our most robust conclusion is 
that planetesimal growth in an initially gravitationally stable turbulent disc is not prohibited 
by turbulence with 10~ 6 < a < 10~ 3 , even in the absence of including dust agglomeration 
enhancement mechanisms. For a > 10~ 2 we do not find much opportunity for growth, even 
for discs with a dead zone, and additional physics to promote dust growth such as large scale 
vortices, or the consequences of mass pile up at R = R tp need to be considered. 

9. Appendix 

Here we derives the various formulas for Tables 3 and 4. 



9.1. Regime 0.1 

In this regime we have H d = H dl = c^fi" 1 from (33) and v c = v C)0 = v t h = \j^ 1 c s from 
(26). Here the dust motion is sufficiently strongly coupled to the bulk motion of the gas that 
the turbulent eddies have only the effect of preventing settling. Accordingly, (3) becomes 



= pE d Q* — = pE d Q x — 9 -p d . (49) 
at y m d V 47T0 d 

Solving (49), the time scale £i(0/) for regime 1 to grow dust grains from (pi to 0/ is: 

= h^^k (*' - **) - l^kfl ^ ^ (50) 

Because the growth rate in this regime falls off strongly with grain size, the time spent in the 
regime depends only on its final size. This allows us to neglect any uncertainty in the initial 
radii of interstellar dust (which we will therefore approximate to be 0). We can also see 
that, as suggested in Sec. 2.1, the actual size dust grains can reach in regime 0.1 is strongly 
limited by time because t ~ 02. 

For the parameters we consider, regime 0.1 will lead into regime 1.1 and thus the former 
ends when t> Cj0 = f c ,i, after which v c = v c> i. Setting v c from (26) and (25) equal we can see 
that this occurs for 

/27m„p2£2£ ml 
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for Epstein drag and for 



<t>f = ^o.i,5 = r£f (52) 



for Stokes drag. 



9.2. Regime 1.1 

When (25) is the dominant value for v c and the dust is too small to settle (H d = H dj i 
from (33)) we are in regime 1.1. In this regime v c is dominated by the smallest scale turbulent 
eddies, approaching the corresponding eddy speed as the dust grows. 

Solving (3) in this regime we find the time scale ti.i(0i,0/) to grow the grains from 0« 
to 0/ is 

t"AM) = ^\fE^(*r) ( 53 ) 



2p V oiVt \ (f) 

for Epstein drag and 



, i± i \ 3^qlmfpPd t-m | 1 1 | ^^q^mfpPd t-ni 1 fr a\ 

*■■'■»<*• *'> = -^ri^i ^i^ta (54) 

for Stokes drag. Note that these values for t\,i are not strongly dependent on 0/, allowing 
us to use the approximations implicit in (25). 

Once the damping time grows long enough such that 

r = ^- (55) 

•Emax 

then the "Goldilocks" turbulent scale exists, v c = i> Cj 2 (24) applies and the regime 1.1 pro- 
gresses to regime 2.1. This occurs for 

h = h,,E EE (56) 



'-'max 



with Epstein drag and for 



3lmfp^gQ>Pdt 



0/ = 01,,S EE ( ) 2 • (57) 

^•^max 



with Stokes drag. 
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9.3. Regime 2.1 

In regime 2.1, the dust grows to couple with progressively larger and larger scale eddies 
that determine the collisional velocity v c . For the range of a and R that we consider, the 
grains still remain small enough to avoid settling. (For smaller a and larger R, even the 
largest eddies would be be too weak to prevent settling and an additional Regime 1.2 would 
occur.) 

In regime 2.1, we have v c ^ from (24) and H^i from (33). Solving (3) we find the growth 
time scale in this regime to be 



for Epstein drag and 



for Stokes drag. 

Regime 2.1 progresses to regime 2.2 when the grains grow large enough that the turbu- 
lence can no longer prevent settling at which point (34) applies. This occurs for 



Accordingly Regime 2.1 ends for 



>2.1,E 



(60) 



(61) 



with Epstein drag and for 



0/ = 02.1,5 = ( 2 J 1 • (62) 



with Stokes drag. 



9.4. Regime 2.2 



Regime 2.2 occurs for dust grains large enough to begin settling out of the gas disc and 
also have their collisional velocity determined by coupling to non-minimum scale turbulence. 
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Here we use v C} 2 from (24) and Hd^ from (34) while approximating 1 — e n T = 1. Solving 
(3) we find the time scale to grow from some 0« to 0/ in this regime to be 

with Epstein drag and 

= 2^ (^) 4 2 - ^ - (^) 4 < M > 

Regime 2.2 ends when the turbulent length scale on which the dust's collisional velocity 
depends reaches the largest turbulent scale. This happens when x — ^m. — XmaX) which 
occurs for 

^nx— (65) 

Accordingly, we find that in the Epstein regime, regime 2.2 ends for 

0/ = 02.2,£ = (66) 



J max 



and 



(67) 



in the Stokes regime. None of (63), (64), (66) or (67) depend on a so the only a dependence 
of the time scales that result from regime 2.2 derives from the initial size condition that 
determines when the regime applies. 



9.5. Regime 3.2 

In regime 3.2 we use v c ^ from (27) while continuing to use H^2 from (34): the dust is 
decoupling from the turbulence while continuing to settle. Writing (3) in this regime we find 

^=p£ d fi 3 T 2 (l-e-w) 2 . (68) 

While the above equation is not simply solvable and in our plots we integrate numerically, 
Regime 3.2 can be divided into two subregimes a and b which admit the approximate forms: 



(69) 
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for r < f2 1 and 

(70) 

for t > n- 1 . 

Using the first time scale approximation we find that the approximate time scales to 
grow from some 0, to 0/ < 0(r = Q^ 1 ) in this regime are less than 

9E 2 

^' wd'-.-) ^ 1 -^ (71) 



for Epstein drag and 



- ( ^-^)~ -r:-^ .- 3 (72) 



for Stokes drag. Using the second time scale approximation, we find that the approximate 
time to grow from some 0j > 0(r = fi _1 ) to 0/ is 

t3 - 2fe - ~ e^" (73) 

Regime 3.2 ends when the dust disc becomes gravitationally unstable, which occurs 
when the Toomre criterion (Binney & Tremaine 1987) 

H, < ^ (74) 
is satisfied. This condition can be rewritten as 

Regime 3.2 then ends when 



(75) 



vrG E 
for Epstein drag and when 

f 3-\/al m fpQc s S g 



3y^C s ^ E s 

0/ = 03.2,1? = ^ \T~ (76) 



for Stokes drag. 
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Table 1. Key Variables 



Variable 


Description 


o 

r 


Dust disc densitv 




Dust grain density 


l J 9 


VJ OjiJ v ill >V LHjllkjll y 




Dnst Hi^r surfflpp Hpnsit,v 


^9 


Gas disc surface densitv 

V_ 1 <_Aji_j LlllJL. tj l_AA X CAi V_; V_X.V_;XXlJX U V 


/ j ,,n 


Gas disc surface density at R = 1AU 


Ta 

' ll 


Grain radius 


r 


Grain size parameter 


77? j 


Grain mass 

VJ A CAiXXX XXXCAilJkJ 


777, „ 

Ilbg 


Mass of a, Eras molecule 


v th 


Dust thermal sneed 


1 1, 


Grain frictional stonninf? time for Dnstein drat? 

V_ 1 A <_Aj A XX XX X V U X V_y XXCAiX l_J U w l_7 lllw U llllv 1 Vl X i l-J l_J U VjI XX vll CAiCL 


TS 


Grain frictional stopping time for Stokes drag 


n 


Keplerian angular velocity 


h 9 


Gas disc scale height 


H d 


Dust disc scale height 


H D 


Dead zone height 


c s 


Gas sound speed 


V 


Sticking parameter 


a 


Shakura-Sunyaev viscosity parameter 


a gg 


Neutral gas collisional cross-section 


CTgd 


Dust-gas collisional cross-section 


&dd 


Grain-grain collisional cross-section 


R 


Distance from the star 


tr 


Total time scale to achieve gravitation instability 



Table 2. Regime Descriptions and Ordering 



Regime Start Condition End Condition Collisional Velocity Dust Height 



Regime 3.2 



Regime 2.2 



Regime 2.1 



T = 



T = 



T = 



T = 



V aQrBcg 



Vacs 

TFT 



Vacs 



n 2 i 



Regime 1.1 Eq. (51), (52) 



T = 



— TC 

(-7T- 



Regime 0.1 Interstellar Dust Eq. (51), (52) 



Note. - The planetesimal growth regimes that occur in our model which progress 
from bottom to top as the dust grains grow. . The collisional velocities are developed in 
equations (26), (25), (24) and (27) and the scale heights in (33) and (34). We use these 
values to calculated growth rates by plugging the collisional and height velocities into eq. 
(3). The beginning and end conditions define the limits of integration for the time scale 
calculations. 
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Table 3. Regime Dust Sizes 



Regime Final 4>e 



Final 4>s 



Regime 3.2b 
Regime 3.2a 
Regime 2.2 
Regime 2.1 
Regime 1.1 
Regime 0.1 



3^/ac s Q S 9 



3S„ 



3S fl 



27m g p 2 d T, 2 g t m n 



( 3^/al mfp nc s S g \ 2 
^ ttG V A Pd) 

( 3lmfp'EgPd \ 2 

V 2 / 

1 

) 



3l m fp'Sgpd \ 2 
2^m 



\ ^Xmax J 



^m g i 2 mf s 2 g np^t m y 



Note. - - Values of bounding the various regimes 
for both Epstein and Stokes drag. Derives from the 
bounding values for r in Table 2. 



Table 4. Regime Time Scales 



Regime Epstein Time Stokes Time 

Regime 3.2b As Epstein 

Regime 3.2a - (<Pf - 

Regime 2.2 ^ (^f {<t>^ ~ <t>?) ^nE § {*7* ~ <t>?) 

Regime 2.1 ^^(4-4) 



2a Y. d tlB 1J "' \ 4>i 

Regime 1.1 ffe^ ($ ) (i " £ 



5 



Regime 0.1 i As Epstein 



Note. - - Equations for the time for growth from = fa to = 0/ in the 
various regimes, presuming p — 1 (optimal sticking), derived by integrating 
Eq. 3. The time scales are inversely dependent on p and are derived in the 
appendix. 
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Fig. 1. — A schematic representation of fractally embedded turbulence. The largest circle is 
an eddy too large to effect the grain, the smallest circles eddies too weak and the medium 
circles the eddy scale with the dominant effect. This is the Goldilocks scale discussed in the 
text. As the grain exits one eddy, there is another for it to enter. 
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R ( AU ) " R (AU) 



Fig. 2. — Size of the dust grains at regime boundaries as a function of the distance from the 
star for pd = lgcm~ 3 . The top pair of panels are for a = 10 -2 while the bottom pair are for 
a = 10~ 4 . The solid curve is the dividing line between Epstein (right) and Stokes (left) drag. 
It lies above the plotted region for the top left panel. The dashed and dot-dashed curves 
represent the grain sizes for which either the collisional velocity of the grains or the dust 
disc scale height changes its dependence on the drag stopping time r (see Table 2). In order 
of increasing grain size, the dashed and dot-dashed curves denote the regime transitions 
between regimes 0.1 and 1.1, 1.1 and 2.1, 2.1 and 2.2, 2.2 and 3.2, and the curve r = 
The highest curve is the grain size at which gravitational instability occurs. 
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R (AU) R (flU) 

Fig. 3. — Total time scales to grow through all regimes from Table 2 up to the gravitionally 
unstable regime for various values of a as a function of distance from the star. The solid 
line is a = 1CT 4 . The short dashed and dot-dashed lines are a of 1CT 2 and 1CT 3 respectively 
while the long dashed and dot-dashed lines are for as of 1CT 5 and 1CT 6 . For any R there is 
a value of a that minimizes the total time scale. We have used pd = lgcm~ 3 and p = 1 for 
these plots. The various time scales vary inversely with p (t ~ 
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Fig. 4. — The time scale for growth through all the regimes to gravitational instability as in 
Fig. (3), but for a disc with live zone a = 1CT 3 and a dead zone with a = 10~ 4 (solid line). 
Also shown are the time scales for discs of uniform a = 1CT 3 (dash-dotted) and a = 1CT 4 
(dashed). The various time scales vary inversely with p (t ~ 
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6000r 




Fig. 5. — The time scales for growth through all the regimes to gravitational instability as 
a function of the distance from the star and for a range of parameters. We use parameters 
from our MMSN (section 2.2) and set the dust grain density pd = lgcm -3 except as listed 
below. The top row halves the dust grain density (p d = 0.5gcm~ 3 ). The middle row doubles 
the surface density of both the gas and dust discs (E ff (lAU) = 2E g0 , = 0.01). The 

bottom row doubles the ratio of dust to gas in the disc (E d /E ff = 0.02). As in Fig. 3 the 
solid lines are a = 10~ 4 , the short dashed and dot-dashed lines are a of 10~ 2 and 10~ 3 
respectively while the long dashed and dot-dashed lines are for a = 10~ 5 and a = 10~ 6 . In 
all panels we have optimal sticking (p = 1). The time scales vary inversely with the sticking 
parameter p (t ~ 1/p). 
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2345678 2345678 
R (AU) R (AU) 



Fig. 6. — Maximum collisional velocity in km/h as a function of distance from the star for 
a live zone a of 10" 3 (left) and 10~ 4 (right). The higher solid lines are the velocities from 
(28) for the live zone while the lower solid lines is the live zone velocity from (41) for a disc 
with a dead zone with Hp = 0.5H g . The long dashed, short dashed and dash-dotted lines 
are the collisional velocities at r = Q^ 1 and the ends of regimes 2.2 and 2.1 respectively. 
Because these curves apply for grains of cm to meter size scales, common experience suggests 
that sold grain destruction from collisions > 30km/h would lead to obliteration rather than 
sticking. This renders further growth impossible for a live zone a > 10~ 3 unless there is a 
contemporaneous dead zone. 
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Fig. 7. — The ratio of the growth time scale %■ to the infall time scale — for the headwind 

O O T VR 

induced infall as calculated in section 6.1.2, as a function of distance from the star. In this 
plot we take p — 1, but multiplying the curve by 1/p characterizes the result for other values 
of the sticking parameter. The discontinuity at R ~ 6AU is due to the change in drag regime 
since Ts ~ 4> 2 whereas te ~ 0. 
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0.01 0.02 0.05 0.1 0.2 0.5 



P 

Fig. 8. — The distance from the star in AU at which the headwind induced infall time of 
Sec. 6.1.2 equals the growth time scale of r, as a function of the sticking parameter p. This 
distance is an approximate measure of the radius inside of which the headwind effect can be 
neglected. 



